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Abstract 

Suppose we are given a set of objects that cover a region and a duration associated with each object. Viewing the 
objects as jobs, can we schedule their beginning times to maximize the length of time that the original region remains 
covered? We call this problem the SENSOR COVER PROBLEM. It arises in the context of covering a region with 
sensors. For example, suppose you wish to monitor activity along a fence (interval) by sensors placed at various fixed 
locations. Each sensor has a range (also an interval) and limited battery life. The problem is then to schedule when 
to turn on the sensors so that the fence is fully monitored for as long as possible. 

This one dimensional problem involves intervals on the real line. Associating a duration to each yields a set of 
rectangles in space and time, each specified by a pair of fixed horizontal endpoints and a height. The objective is 
to assign a bottom position to each rectangle (by moving them up or down) so as to maximize the height at which 
the spanning interval is fully covered. We call this one dimensional problem RESTRICTED STRIP COVERING. If we 
replace the covering constraint by a packing constraint (rectangles may not overlap, and the goal is to minimize the 
highest point covered), then the problem is identical to DYNAMIC STORAGE ALLOCATION, a well-studied scheduling 
problem, which is in turn a restricted case of the well known problem STRIP PACKING. 

We present a collection of algorithms for RESTRICTED STRIP COVERING. We show that the problem is NP-hard 
and present an 0(log log n)-approximation algorithm. We also present better approximation or exact algorithms for 
some special cases. For the general SENSOR COVER PROBLEM, we distinguish between cases in which elements 
have uniform or variable durations. The results depend on the structure of the region to be covered: We give a 
polynomial-time, exact algorithm for the uniform-duration case of RESTRICTED STRIP COVERING but prove that 
the uniform-duration case for higher-dimensional regions is NP-hard. Finally, we consider regions that are arbitrary 
sets, and we present an 0(log n) -approximation algorithm for the most general case. 

1 Introduction 

Sensors are small, low-cost devices that can be placed in a region to monitor local conditions. Distributed sensor 
networks have become increasingly more popular as advances in MEMS and fabrication allow for such systems that 
can perform sensing and communication. How sensors communicate is a well-studied problem. Our main interest 
is: Once a sensor network has been established, how can we maximize the lifetime of the network? It is clear that 
the limited battery capacities of sensors is a key constraint in maximizing the lifetime of a network. Additionally, 
research shows that partitioning the sensors into covers and iterating through them in a round-robin fashion increases 
the lifetime of the network [1, 3, 8, 9]. 

Definitions. Let S = {si, . . . , s n } be a set of sensors. Each sensor s can be viewed as a point in some space with 
an associated region R(s) of coverage. For every point x <G R(s), s is said to be live at x. Let U be the region to 
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Table 1: Summary of results. 
Uniform duration 



Variable duration 



Intervals 



exact in P 



NP-hard, O (log log n) -approx. 



Rectangles, Disks, . . . 



NP-hard, 0(ln(n/L))-approx. 



NP-hard, 0(logn)-approx. 



Arbitrary sets 



logn-hard to approx., 
0(logn)-approx. [5] 



log n-hard to approx., 
0(log n) -approx. 



be covered by the sensors. U is covered by a collection 1Z C S of sensors if U C [J sen R(s). We often refer to 
TZ as a feasible cover. Every sensor s £ S can be active for a finite duration d(s). Let d m [ n = min se 5 d(s), and 
tax = max se5 d(s). 

Problem (Sensor Cover). Compute a schedule S of maximum duration T, in which each sensor s £ S is assigned 
a start time t{s) > 0, such that any x £ U is covered by some active sensor at all times < t < T. That is, for all 
x £ U and < t < T, there is some s € S with x £ R(s) and t(s) <t< t{s) + d{s). 

A sensor is redundant in a schedule S if it can be removed without decreasing the duration of S. A schedule 
with no redundant sensors is a minimal schedule. By removing the redundant sensors, any optimal schedule can be 
converted to a minimal schedule of the same total duration. Therefore it suffices to consider only minimal schedules, 
which may not utilize all sensors. As a convention, we set t(s) = oo if s is unused. 

Prior work on the SENSOR Cover problem has focused solely on the case where the regions R(s) are arbitrary 
subsets of U and the durations are all identical. This assumption yields a packing constraint, and the problem reduces 
to partitioning the set of sensors into a maximum number of valid covers. This problem is known as Set Cover 
Packing and is Inn-hard to approximate, with a matching upper bound [5]. 

In practice though, these assumptions appear overly constraining. Sensors will have arbitrary durations and typ- 
ically define geometric regions of coverage: intervals, rectangles, disks, etc. In this paper, we will consider these 
classes of problems. In the Restricted Strip Cover problem, all the R(s)'s are intervals in one dimension. In this 
case the problem is equivalent to sliding axis-parallel rectangles vertically to cover a rectangular region of maximum 
height. Thus, Restricted Strip Cover resembles the Dynamic Storage Allocation (DSA) problem [6], 
albeit in a dual-like fashion. In the Cube Cover problem, the i?(s)'s are axis-parallel rectangles, and the problem is 
akin to sliding cubes vertically in the z-dimension. We also consider SENSOR COVER, when the R(s)'s are arbitrary 
subsets of a finite set U of size \U\ = 0(n), with varying durations (in contrast to Set COVER PACKING). 

In general, a schedule may activate and deactivate a sensor more than once. We call this a preemptive schedule. A 
non-preemptive schedule is a schedule in which each sensor is activated at most once. In this paper we only consider 
the non-preemptive problem. We have some preliminary results for the preemptive case, but more research is needed 
to gain a better understanding of the differences. 

Our results. We show that most variants of SENSOR COVER are NP-hard, and we study approximation algorithms. 
For any point x £ U, let L(x) = ^2 seS s n ve3tx d(s) be the load at x . Define the overall load L = mm x L(x). 
We write Lx (rsp., Lx(x)) for the load of any subset X of sensors (rsp., at x). Letting OPT denote the duration of 
an optimal schedule, a trivial upper bound is OPT < L. All our approximation ratios are with respect to L. That 
OPT < L allows the assumption that d max < L, because durations exceeding L contribute to neither load nor OPT. 

Table 1 summarizes our results. The most interesting case is the Restricted Strip Cover problem, for which 
we give an 0(log log ^-approximation algorithm. After reviewing some related work, we discuss RESTRICTED 
Strip Cover in Section 2, Cube Cover in Section 3, and the general Sensor Cover problem in Section 4. 

Related Work. Set Cover Packing was studied by Feige et al. [5]. They considered the Domatic Number 
problem, where the goal is to maximize the number of disjoint dominating sets on the set of vertices of a graph. A 
dominating set in a graph G = (V, E) is a set V C V of vertices such that every v £ V is either contained in V 
or has a neighbor in V'. Feige et al. show for every e > 0, the DOMATIC NUMBER problem is hard to approximate 
within a factor of (1 - e) ln|V|. In the proof of their hardness result, Feige et al. use the Set Cover Packing 
problem, whose goal is to maximize the number of disjoint set covers, given a set of subsets S of a base set U. Note 
that the Set Cover Packing problem is a combinatorial version of our problem, with each subset being a region, 
and each (sensor) region having unit duration. Feige et al. also give a randomized Inn-approximation algorithm, 
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which they derandomize. Another key feature of this work is that it showed the first maximization problem proved to 
be approximable within polylogarithmic factors but no better. 

The practical motivations for studying this problem have led to the development of numerous heuristics. Slijepce- 
vic and Potkonjak [9] introduce the Set K-Cover problem, where they are given a set of subsets of a base set and 
an integer k and ask if it is possible to construct at least k disjoint set covers. They re-prove an NP-hardness result for 
this problem, probably unaware of the Feige et al. result. They also present a heuristic that selects mutually exclusive 
sets of nodes, where each set completely covers the desired region. 

Perillo and Heinzelman [8] study a variation of this problem, where they want to maximize the lifetime of a multi- 
mode sensor network. They compute all possible feasible covers and then translate their problem instance into a graph. 
Each sensor and feasible cover becomes a node. Sensors are connected to a feasible cover if they are contained in that 
feasible cover. They use linear programming to model additional energy constraints and solve the maximum flow 
problem on this graph. Although they solve the problem optimally, their solution can be exponential in the problem 
instance. Dasika et al. [3] also compute all possible feasible covers and develop heuristics for switching between these 
covers in order to maximize the lifetime of their sensor network. 

Abrams et al. [1] study a variation of the problem where they are given a collection of subsets of a base set and a 
positive integer k > 2. Their goal is to partition the subsets into k covers, where the area of coverage, defined as the 
cardinality of a set, is maximized across all k covers. They give three approximation algorithms for this problem: a 
randomized algorithm, a distributed greedy algorithm, and a centralized greedy algorithm. Their randomized algorithm 
partitions sensors within 1 — - of the optimal solution. Their distributed greedy algorithm gives a | -approximation 
ratio. Their centralized greedy algorithm achieves an approximation factor of 1 — K They also prove a y|-hardness 
result for their problem. 

To our knowledge, the Restricted Strip Cover problem has not previously been considered. Some of the 
closely related problems are well studied, however. If we replace the covering constraint by a packing constraint 
(rectangles may not overlap, and the goal is to minimize the height of the highest point covered), then the problem is 
the same as DYNAMIC STORAGE ALLOCATION [6, Problem SR2], for which there is a (2 + e)-approximation [2]. 
If we further allow rectangles to move both vertically and horizontally, then the problem becomes STRIP PACKING, 
which has a (1 + e)-approximation (up to an additive term) [7]. 

2 Restricted Strip Cover 

Consider an instance S of Restricted Strip Cover (RSC). For ease of presentation, we define R(s) as a semi- 
closed interval [£(s),r(s)) for each s £ S. We assume without loss of generality that all interval coordinates are 
non-negative integers in the range [0, 2n — 1] and U = Uses ^( s ) = [0; 2n — 1), because there are at most 2n 
distinct interval endpoints of sensors. It is convenient to view scheduled sensors as semi-closed rectangles in the 
plane, with their intervals along the x-axis and durations along the y-axis. Thus a valid schedule S of duration T is 
one in which any point (x,y) in the sub-plane U x [0,T) is covered by some sensor s; i.e., l(s) < x < r(s) and 
t( s ) < V < t( s ) + d(s). The problem is equivalent to sliding axis-parallel rectangles vertically to cover a rectangular 
region of maximum height. Therefore, in this section we use the terms "sensor" and "rectangle" interchangeably. We 
say two or more rectangles overlap if they cover some common point. When discussing multiple schedules, we write 
ts(s) to denote the start time of s in some schedule S; otherwise we omit the subscript. 

We also assume that all durations are positive integers. We define level j of any schedule to be the horizontal slice 
of sensors that cover points at y-coordinates in [j — Let S be some schedule of S. A gap is a point p such that 
no sensor covers p. For any i € U, define M(S, i) to be the greatest y-coordinate j such that no gap exists below j at 
i; i.e., M(S, i) = max{j : Vf < j, 3s € S, s covers (i, j')}. Then the duration of S is M(S) = mini M(S, i). 

Our main result for RSC is an O (log log n) -approximation. En-route we obtain better results for some special 
cases: (i) a simple, exact algorithm if all sensors have the same duration (Section 2.1); (ii) an exact, dynamic pro- 
gramming algorithm, which runs in poly(n) time if L = 0(logn/ loglogn) (Section 2.3); and (iii) 1/(1 + e)- 
approximations if L = 0(d min log nj loglogn) or L = £!(d raax logn • (1/e) 4 min{l/e, log(d raax /d min )}) (Sec- 
tion 2.4). 

Our main technique builds on the method outlined by Buchsbaum et al. [2] in their algorithm for DYNAMIC 
Storage Allocation. Their approach does not apply directly, because the covering constraint poses different 
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Figure 1: (From Buchsbaum et al. [2].) A set of sensors (in the form (£(•), r(-), A = (0, 1, 3), B = (0, 3, 1), 

C = (1, 2, 2), D = (1, 4, 1), E = (2, 3, 1), F = (2, 5, 1), G = (3, 4, 2), and H = (4, 5, 3). The shaded region is a 
gap. In this example, L = 4 but OPT = 3, which can be realized by sliding G down so that t(G) = 1. 

challenges to solving RSC than does the packing constraint to DYNAMIC STORAGE ALLOCATION. We will adapt 
parts of their method to obtain our bounds. 

2.1 Uniform-Duration Sensors 

If all sensors have the same duration, a simple greedy algorithm gives an exact solution of duration L. Define Si = 
{s G S : s is live at i}. Assume by scaling that all sensors have unit duration. We proceed left-to-right, starting at 
i = and constructing a schedule S while maintaining the following invariants after scheduling sensors in Sf. (i) no 
sensors overlap at any ^-coordinate > i, and (ii) M(S, i) = L. 

When i = 0, select any L sensors that are live at 0, and schedule them without overlap, establishing the initial 
invariants. Assuming the invariants are true at i, schedule Si+i as follows. If there are no gaps at i + 1, we are done, 
as the invariants extend to i + 1. Otherwise, assume there are k > unit-duration gaps at i + 1. The invariants imply 
that at least k sensors in Si+i remain unscheduled, which can be used to fill the gaps while preserving the invariants. 

2.2 Hardness Results 

When sensors have variable durations, the problem becomes NP-hard. To prove this, we exploit an identity to DY- 
NAMIC Storage Allocation in a special case. An instance of DS A is like one of RSC, except that the instance 
load is defined as the maximum load at any x-coordinate, and the goal is to schedule the sensors (jobs, in DYNAMIC 
Storage Allocation parlance) without overlap so as to minimize the makespan. If the load is equal for all x- 
coordinates, then OPT = L for either problem implies a schedule that is a solid rectangle of height OPT = L. For 
RSC, it further implies that the schedule is non-overlapping; for DS A, this implication is redundant. 

Stockmeyer proves that determining if there exists a solution to DS A of a given makespan is NP-complete [6, 
Problem SR2]. His proof reduces an instance of 3-PARTITION to an instance of DSA with uniform load and, in 
fact, all durations in {1, 2}, such that OPT = L if and only if there is a solution to the 3 -PARTITION instance. (See 
Appendix A for details.) Thus he proves that given a DSA instance with uniform load, determining if OPT = L is 
NP-complete, even if durations are restricted to the set {1, 2}, and by the above identity, the same is true for RSC. 

To establish a gap between OPT and L, consider the example in Figure 1, in which L = 4 but OPT = 3. Scaling 
the durations shows that no approximation algorithm can guarantee a ratio of better than 4/3 with respect to L. 

2.3 A Dynamic Programming Solution for Small L 

We give a dynamic program to answer the question: Is there a schedule £ such that M (S) = T for a fixed T? In the 
following, we ignore portions of sensors that extend above level T in any schedule. 

Define 6><j = Uo<fe<i^ fe - Consider some schedules of <Sj_i and Si of Si such that M(Si-±,i — 1) = 
Al (Si, i) — T. We say that and Si are compatible if (i) for all s £ Si-i n Si, ts i _ 1 (s) = ts; (s); and (ii) for all 
j € [0, T), (i, j) is covered by Si-i or S^ The first condition stipulates that any sensor in both schedules must have the 



4 



same start time in each; the second requires a sensor in Si to be scheduled to cover each level at which coverage stops 
at i — 1 in Si-i. For each i, we populate an array d indexed by possible schedules of Si. For any Si, define Ci[Si] = 1 
if there is a schedule S of S<i that respects Si and has M(S, x) =T for < x < i; and Ci[Si] = otherwise. Then 
Cj[Si] = 1 if and only if M(Si, i) — T and there exists some schedule Si-i of <Sj_i such that Cj_i [S^-i] = 1 and 
Si-i is compatible with Si. For i = 0, Co [So] = 1 f° r precisely those schedules So of So that have M(Sq, 0) = T. 
The dynamic program then populates the arrays Cj in increasing order of i, by checking all schedules of Si for each i. 
Ultimately we check if there is some schedule S^n-i of S2 n -i sucn that C2n-i[S2n-i] = !■ 

First Analysis. For a schedule Si of Si, consider the union of the rectangles of Si, and denote by d (Si) the vertical 
boundaries of this union. If Si is part of a minimal schedule S of duration T, then any rectangle of Si must cover some 
point on d(Si) that is covered by no other rectangles in Si. Thus |Sj| < 2T, because d(Si) has total length 2T. 

Now we can analyze the dynamic program, which we restrict to consider only minimal schedules. The number of 
sensors in any schedule of Si is at most 2T, so there are at most {^)T 2T possible schedules of Si, as each potential 
set of r\ sensors can be scheduled in T v ways. Each schedule of Si must be checked for compatibility against each 
schedule of <Si_i, and checking compatibility of a pair of schedules takes 0(T) time. Hence the time to run the 
whole dynamic program is 2n((™ T )T 2T ) 0(T) = (nT)°^ = (nL)°( L l To determine OPT, we run the dynamic 
program for each of the L possible values of T, which does not affect the overall asymptotics. 

Partitioning the Dynamic Program. Now we restrict the x-coordinates on which we have to run the dynamic program 
to those with relatively few live sensors. Let X = {i : |<Sj| < 5T}. We claim that S has a schedule of duration T if 
and only if S has a schedule S such that M(S, i) >T for any i £ X. We prove the "if" part; the "only if" part is clear. 

Assume that there is a minimal schedule S of duration T that only covers X. We show how to schedule the sensors 
not used in S to cover all ^-coordinates. Consider any maximal interval X of x-coordinates not in X. At most AT 
sensors from S are live at any i e X, because any such sensor is also live at either min(X) — 1 or ma,x(X) + 1, and 
at most 2T are live at either one. By construction, there are at least 5T sensors live at any i € X, so there are at least 
5T — AT = T sensors live at i that are not used by S and hence are available, which suffice to cover all the levels at i. 
If such a sensor s should also be live at another i' € X (or another i' in another X'), it reduces by one both the number 
of potential uncovered levels and the number of available sensors live at i' , so enough sensors will remain at i'. 

Therefore we need only run the dynamic program on the cc-coordinates in X. This takes only 2n ■ T°^ time, 
because there are fewer than 5T sensors live at any i € X . Thus we have: 

Theorem 2.1. RSC can be solved exactly in time 2n ■ L°( L \ 

Corollary 2.2. RSC can be solved exactly in poly(n) timeifL < c-log n / log log n for some constant c small enough. 

Using a standard trick, a PTAS follows directly by truncating durations appropriately. 
Corollary 2.3. There is a PTAS for RSC if L < c ■ <i m i n log nj log log n for some sufficiently small constant c. 

2.4 Approximation Algorithms via Grouping 

In this section, we give approximation algorithms via the grouping technique, which is similar to the boxing technique 
of Buchsbaum et al. [2]. We know that the load L is a natural upper bound on OPT, and L = OPT when all sensors 
have the same duration. The basic idea of grouping is to group shorter sensors into longer, virtual sensors until all the 
sensors have equal duration, at which point the greedy algorithm is invoked. Essentially we must ensure that the load 
does not decrease too much during the process, which is the central component of our algorithms. 

Grouping Sensors. A grouping of a set Y of sensors into a set of groups G is a partition of Y into \G\ subsets, each 
of which is then replaced by a rectangle that can be covered by the sensors in the group. The duration of a group 
is defined to be the duration of the rectangle that replaces it. That is, these rectangles can be viewed as sensors in a 
modified instance. Then Lq (rsp., Lq(i)) is defined to be the load of the groups (rsp., at i). Note that Lq(i) < Ly(i), 
since portions of the sensors in a group that are overlapped or outside the rectangle are not counted in Lq(i). In the 
following, we give procedures to group a set Y of sensors of unit duration into G such that Lq(i) is not much smaller 
than Ly(i) for any i. All of the grouping procedures in this section run in polynomial time. 
First, we give a grouping of a set of sensors that are all live at a fixed x-coordinate. 



5 




Figure 2: (From Buchsbaum et al. [2].) (a) Four sensors, (b) The sensors of (a) viewed as (x, y) points, (c) Grouping a 
set of sensors with D = 2 and e = 1/2. The rectangle Ri contains the set Y of sensors. The sensors are first partitioned 
into alternating vertical and horizontal strips of D\l/e] = 4 each. Within each strip, the sensors are grouped (dotted 
lines) into groups of D = 2. The groups that intersect the line x = i are shaded. 



Lemma 2.4. Given a set Y of unit-duration sensors, all live at some fixed x-coordinate x$, an integer group-duration 
parameter D, and a sufficiently small positive e, there is a set G of groups, each of duration D, such that for any i, 

L G (i) > L Y (i)/{l + e) - 4D\l/e\. 

Proof. It is convenient to view a sensor s as a point (£(s), r(s)) in the plane. Note that all sensors live at xq are inside 
the rectangle R X(I = {(x,y) : x < xq < y} (Figures 2(a)-(b)). First we partition the sensors of Y into strips by 
repeating the following as long as sensors remain. 

(1) Create a vertical strip containing the at most D\l/e] sensors that remain with the smallest £(■) values. 

(2) Create a horizontal strip containing the at most D \1 /e] sensors that remain with the largest r(-) values. 
Now for every vertical strip of Y, take the sensors in order of decreasing r(-) value in groups of size D (we 

may discard the last < D sensors in the last strip). Similarly, for every horizontal strip, take the sensors in order of 
increasing £(■) value in groups of size D (we may discard the last < D sensors in the last strip). Replace each group 
X with a larger rectangle sx with £{sx) = max se x £{s), r(sx) = min se x r ( s ), an d d(sx) = J2sex d( s ) = D. 

Consider any i < xo (the case i > xo is symmetric), and examine Figure 2(c). All sensors live at i are inside 
the rectangle Ri = {(x,y) : x < i < y}. Assume that the line x = i intersects k horizontal strips; then Ri entirely 
contains at least k — 1 vertical strips, so Ly{i) > (k — l)D\l/e] . For any group completely inside Ri, it contributes 
D to both Ly(i) and La(i); for any group completely outside Ri, it does not contribute anything to either Ly(i) or 
Lg{i). So only the groups in the k horizontal strips and the single vertical strip intersected by the line x = i contribute 
to the difference, that is, Ly(i) — Lg{i) < kD + D\l/e] + D, where the last term accounts for the fewer than D 
sensors that we did not group in the last strip. Therefore, 

L>G(i) > L Y (i) — (k — 1)D — (2 + \l/e\)D > (1 — e)L Y (i) — 2D\l/e\ > ^-L-L y (i) - 2D\l/e], 

for any e < 1/2. Replacing e with e/2 gives the desired result. □ 

Next we use Lemma 2.4 to group all sensors of unit duration. 

Lemma 2.5. Given a set Z of unit-duration sensors, an integer group-duration parameter D, and a sufficiently small 
positive e, there is a set G of groups, each of duration D, such that at any x-coordinate i, 

L G (i) > L z (i)/(1 + e) - 0(logn • D/e). 
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Figure 3: A bad example 

Proof. Build an interval tree T on the ^-projections of the rectangles of Z. For each node v of T, let Z v be the set of 
sensors associated with v. All the sensors of Z v are live at a fixed x-coordinate, namely the dividing line at v, and thus 
we can apply Lemma 2.4 to Z v for each v. 

Consider any x-coordinate i. For any two different nodes u, v of T at the same level, the sensors in Z u and the 
sensors in Z v do not overlap. Thus the sensors live at i are distributed to at most O(logn) nodes in T, because 
the interval tree has height 0(log?i). By Lemma 2.4, we have La(i) > Y^ V £t(Lz v (*)/(! + e ) — 4-D|~l/e~|) = 
L z (i)/{l + e)-0(\ogn-D/e). □ 

Remark. The 0(log n) factor in the error term of Lemma 2.5 cannot be removed. Consider grouping the example in 
Figure 3 with D = 2. First, at least half of sensor A has to be wasted, because it is either grouped with some sensor in 
the left half or some sensor in the right half. Assume it is grouped with some sensor in the left half. Then by a similar 
argument, sensor B is cut in half, and one of the halves has to be wasted. Ultimately, we can find an x-coordinate i 
without a single group covering it; i.e., Lc(i) = 0, but Lz(i) = fi(logn). 

The Algorithm. Let e be a sufficiently small error parameter, and let D = d max \l/e\. 



Algorithm 1 Approximation algorithm via grouping 

(1) Truncate each sensor of duration d to [(1 + e) ], where (1 + e) k < d < (1 + e) k+1 for some integer k. Let X 
be the set of truncated sensors. 

(2) For each d = \(1 + e) fe ] , k = [log 1+e <i m i n J , . . . , [log 1+e d max — 1], do the following. Let Xd denote the 
set of truncated sensors of duration d. Scale each sensor in Xd down by a factor of d, apply Lemma 2.5 with 
group-duration parameter \D/d~\ and the given e, and then scale the obtained groups back up by d. 

(3) Let G be the set of rectangles obtained from Step (2). Truncate them so that they all have duration exactly D. 
Call the resulting set of rectangles G . 

(4) Apply the greedy algorithm to G . 



Theorem 2.6. For any sufficiently small positive e, Algorithm 1 runs in poly(n, 1/e) time and gives a schedule of the 
RSC problem with duration at least L/(l + e) — O (d ma x logn -1/e 3 log(d max /(i m in)). 

Proof. We will show that the truncating and grouping do not decrease the load at any i excessively. 

By Lemma 2.5, Step 2 produces a grouping Gd of Xd of duration \D/d~\ d such that at any i, Lc d {i) > Lx d (i)/(l+ 
e) — 0(logn • D/e). Summing over all d, we have 

T 1 ... /Z?l0gnl0g(d m ax/dmin)\ 1 T , „ / dmax log(d max /d min ) log 11 

Law > ——L x (i) - O — = —— L x {i) - O ; 

1 + e \ elog(l + e) J 1 + e \ e A 

Truncating the sensors in Step (1) decreases their durations by at most 1 + e, so Lx (i) > Truncating the 

groups in Step (3) decreases their durations by a factor of at most 1"-° A'l d < D±d <- ^ _|_ ^ t00 j; mce > yq^, 

we have Lqi(i) > L(i)/(1 + 7e) — O (d max logn • 1/e 3 log(<i max /<i m i n )). Finally, applying the greedy algorithm in 
Step (4) yields a schedule of duration mirii Lc{i) > L/(l + 7e) — 0(d max logn- 1/e 3 log(d max /d m i n )). Replacing 
e with e/7 gives the desired result. □ 
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Corollary 2.7. There is a constant c such that for any small enough positive e, the algorithm gives a schedule of 
duration at least L/(l + e) for any L > c ■ <i max logn • 1/e 4 log(d max /d m i n ). 

An Alternative Algorithm. By bootstrapping Steps (l)-(3) of Algorithm 1, we can replace the 0(log(d ma x/dmin)) 
factor with 0(1 /e), leading to the following result. 

Theorem 2.8. For any sufficiently small positive e, there is an algorithm that runs in poly(n, 1/e) time and gives a 
schedule to the RSC problem with duration at least L/(l + e) — O (d m ax logn • 1/e 4 ). 

Proof. We are going to apply Steps (l)-(3) of Algorithm 1 repeatedly, grouping the smaller sensors so as to increase 
d m in until log(d ma x/rfmin) becomes small enough that we can apply Theorem 2.6 to the resulted rectangles. 

For ease of presentation, we assume that 1/e is an integer. Let r denote the ratio d max /dmin- Assume first that 
log r > 1/e, and set p = e/ log r and D = |~ n 4 <i m ax] ■ Apply Steps ( 1 )-(3) of Algorithm 1 to S s , the set of sensors of 
duration at most d' max = \pD] , with group duration D and error parameter p. This yields a set of rectangles G s of 
duration D such that for any i, 

1 + p V V 

> — L s ,(i) — 0(p 2 d max logn\og(p 5 r)) > — — — L s ,(i) - ClC ^maxlog"- ^ ^ . g SQme constant ) 
1 + p 1 + p log r 

Now consider G s as a set of sensors and the new problem instance S 1 = G s U (S \ S s ). Its load at i is 

T ,. w L(i) cie 2 d max logn 

As'W > T~i ] ■ 

1 + p log r 

Moreover, the new minimum duration of this problem instance is at least d' max , and the maximum duration remains 
d max , so the new ratio is r' < < = lo ^ B r < log 10 r, since log r > 1/e. For e sufficiently small, we have 

r' < \[r\ hence log r' < \ log r. 

Next repeat the procedure above, each time using new error parameter p! = e/logr', until it yields a problem 
instance S* with minimum duration <i* lin for which r* = d max /d* llin is such that logr* < 1/e. Let ro, . . . , = r* 
be the sequence of ratios and Lo(i) = L(i), L\(i), . . . , Ljt(z) = L*(i) be the sequence of loads. We have 

r */.v . 1 t c\ ci£ 2 dmaxlogn 

l + e/logrfe_i logrfc_i 

1 / 1 f c 1 e 2 d meix \ogn\ cie 2 d max logn 



1 + e/logrfc_i \1 + e/logr fc _ 2 logr fe _ 2 / logr- fe _ 



a=0 v 7 / i=0 

1 2 

> ; L(i) cie 2 <i max logn (c 2 is some constant.) 

l + C2e/logr* logr* 

> L{i)/{l + 2c 2 e 2 ) -4cie 3 d max logn. 

Let L* = min p Lg* (i). Finally, apply Theorem 2.6 to S*, which yields a schedule of duration at least 

1 _J* _q( rf max log r* log n \ > 1 £ _ q / rfmaxlogn 



1 + e V e 3 J ~ 1 + c 4 e 

for some constant C4. Replacing e with e/c4 gives the desired result. □ 

Corollary 2.9. There is a constant c, such that for any small enough positive real e, the algorithm gives a schedule of 
duration at least L/(l + e)for any L > c ■ d max logn • 1/e 5 . 



2.5 An O (log log n) -Approximation for Arbitrary L 

Theorem 2.8 yields a good approximation only when rf max is small. To extend this, we separate tall rectangles from 
the short ones and handle the former individually. Henceforth, we will analyze the approximation ratio asymptotically, 
and we assume that all durations are powers of 2, which at worst halves the duration of the schedule. 

Let d max = 2 P , and £ = \2 log log n] . We partition S into I + 1 subsets TZq , . . . , TZe . The first subset TZq consists 
of all sensors of duration at most 2 p ~ e ; then for each p — £ + 1 < k < p, we put all sensors of duration 2 k into one 
subset. We call TZq the small subset and the rest large subsets. For any x-coordinate i, we compute Ls in ji k (i), the 
load of Si n IZk at i, for k = 0, . . . , I. Let m(i) = argmax fe Lg ir[ fi k (i); i.e., the sensors in 5, (~1 TZk have maximum 
load when k = m(i). Break ties arbitrarily. It is easy to see that the load of <Sj ft R m (i) is at least Lj {I + 1). 

Next, for each k, we use TZk to cover all the ^-coordinates i where m(i) — k, for a duration of fl(L/ log log n). 
For a large subset TZk, k > 1, because all its sensors have the same duration, we can use the greedy algorithm to find a 
schedule of duration at least L/(£ + 1) = il(L/ log log n). For the small subset TZq, we use Theorem 2.8 with e = 1. 
Since the sensors in TZq have maximum duration 2 p ~ e < L/2 e , Theorem 2.8 yields a schedule of duration at least 

2(£+l) \2* \log\ognJ \log 2 n & J \log\ognJ 

Theorem 2.10. There exists a polynomial-time 0(log log n)-approximation algorithm for the RSC problem. 

3 Cube Cover 
3.1 Hardness Results 

When the i?(-)' s ar e axis-aligned rectangles and U is a two-dimensional region, the problem is NP-hard even when 
the sensors have uniform duration, in contrast to the uniform-duration case for Restricted Strip Cover. We use 
a reduction from an instance of NAE-3SAT with n variables and m clauses to an instance of Cube Cover. 

An instance / of NAE-3SAT is a set U = {iti, 112, u n } of variables and a collection C = {C\, C2, C m } of 
clauses over U, such that each clause C; € C has |Cj| = 3. The problem is to determine if there a truth assignment for 
U such that each clause in C has at least one true literal and at least one false literal [6]. A key property of NAE-3S AT 
is that if X is a satisfying assignment for an instance / of NAE-3SAT, X is also a satisfying assignment of /. 

Given /, we construct an associated graph G(I), with vertices for each variable and each clause. We draw and edge 
between a clause vertex and a variable vertex if the variable appears in the clause. The graph is drawn on a planar grid 
within a bounding box U. From G(I), we construct an instance S(I) of Cube Cover that has a schedule of duration 
2 if and only if I is satisfiable. If / is unsatisfiable, S(I) has a schedule of duration 1. We describe the construction in 
more detail below: 

Variables. Each variable will be represented by a collection of rectangles that cover a square grid. Each rectangle 
covering the variable gadget has unit duration. The rectangles that cover the variable gadget are shown in 
Figure 4, arranged in the "true" and "false" encodings. Using a mixture of the rectangles in the "true" and 
"false" encodings leads to a suboptimal schedule; such a mixture is called an improper cover. 

Pipes. We connect each variable to each clause that contains it via pipes. The pipes are drawn using polygonal lines 
in the plane. A pipe corresponding to a positive (rsp., negative) occurrence of a variable in a clause leaves via 
the left (rsp., right) side of the variable gadget. All rectangles that cover a pipe have unit duration. A pipe, as it 
is drawn within the bounding box, is the leftmost image in Figure 5. To the right of the pipe are the two proper 
configurations for covering the pipe. The first of the configurations shows how the pipe will be covered when the 
clause to which it connects is satisfied by the variable assignment. We refer to this as an "on" signal. The second 
configuration shows how the pipe is covered otherwise. We refer to this as an "off" signal. Using a mixture of 
rectangles in the two configurations (or an improper cover) will lead to a suboptimal schedule. Figure 6 shows 
how the variables and pipes connect. 
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Figure 4: Variables Figure 5: Pipes Figure 6: Variables and Pipes 

Clauses. A clause is represented by a square grid at the end of a pipe. A clause is covered if a variable contained 
in it is satisfied. All other parts of the bounding box are covered with two of the same rectangles, each with 
unit duration so we can guarantee that every point in U, not part of the variable, pipe, or clause gadgets, can be 
covered for a total of 2 time units. 

In the worst case, each of the variables require a square grid of size 0(m), hence the entire grid must contain 
0(nm) rows and 0(nm) columns. Our construction requires 0(n 2 m 2 ) unit squares and thus can be performed in 
polynomial time and space in the size of the input. 

Proposition 3.1. A variable gadget can be covered for two time units using proper covers. If an improper cover is 
used for d time units, then the gadget can only be covered for 2 — d time units. 

Proof. Each point in the variable gadget is covered by exactly two rectangles. If proper covers are used, the bounding 
rectangle can be covered for two units (one for each cover). Note that this is the only way to partition the rectangles 
into disjoint feasible covers. 

Suppose an improper cover is used for d units and the two proper covers are used for x, y time units. Thus the total 
time for which this rectangle is covered is x + y + d. Since at least one rectangle from each proper cover must be in 
an improper cover, x + d < 1 and y + d < 1, and therefore x + y + d<2 — d. □ 

Proposition 3.2. The clause square is covered if a variable contained in it is satisfied. 

Proof. If clause C contains variable u and variable u is satisfied in that clause, this means the configuration shown in 
the middle image of Figure 5 is chosen so the clause grid is covered. □ 

Proposition 3.3. If X is a satisfying assignment for NAE-3S AT, then X is a satisfying assignment as well. 

Proof. If X is a satisfying assignment for NAE-3SAT, then every clause contains at least one variable, u, that is true 
and one variable, v, that is false. The complement, X, makes u false and v true. X ensures that every clause contains 
at least one true variable and one false variable and thus is a satisfying assignment. □ 

Lemma 3.4. If I is satisfiable, then the instance o/Cube Cover has a schedule of duration 2. 

Proof. Let X be a satisfying assignment of /. For each variable that is set to 1, choose the "1" orientation for the 
corresponding variable gadget and the corresponding "on" signal in the pipes connecting this variable to the clauses 
that contain a positive occurrence of the variable. Choose the "off" signal in pipes connecting this variable to the 
clauses that contain a negative occurrence of the variable. We repeat this process for variables set to "0". 

Since X is a satisfying assignment, each clause grid will be covered. Moreover, since X is also a satisfying 
assignment, we can repeat this process for another time unit. Note that each pipe is used exactly once to send a "on" 
signal and once to send a "off" signal, so it can be covered for two time units. Likewise, each variable is used exactly 
once as "1" and once as "0", so it can also be covered for two time units. □ 

Lemma 3.5. If I is unsatisfiable, then the instance o/CUBE Cover has a schedule of duration at most 1. 
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Proof. If all variables are covered using proper covers, we obtain a valid assignment, so some clause rectangle must 
remain uncovered (since I is unsatisfiable). In order to cover this clause rectangle, we must use an improper cover. 
Our construction requires that an improper cover creates an overlap either in the variable gadget or the pipes. Since 
the load at every point in the variable gadget and the pipes is 2 and a rectangle must be used in its entirety or not at all, 
the maximum schedule has duration 1 . □ 

Theorem 3.6. Cube Cover is NP-hard and does not admit a PTAS, even with uniform duration. 

Proof. NP-hardness follows from Lemma 3.4 and Lemma 3.5. 

Assume we have a PTAS V for CUBE COVER. On input (<S(I), e), where e > and S(I) is an instance of CUBE 
COVER induced by the above construction, V would output a solution with duration T, where T > (1 — e) ■ OPT. 
Setting e = 0.25, T > 1.5 when OPT = 2, and 0.75 < T < 1 when OPT = 1. Thus we can use V to distinguish 
between satisfiable and unsatisfiable instances of NAE-3SAT. □ 



3.2 Rectangles With Unit Duration 

We consider approximation algorithms for Cube Cover if all sensors have unit duration. First we prove a technical 
lemma, which actually holds for arbitrary sets. 

Lemma 3.7. Let U be finite set with m elements. For each s£5, R(s) is an arbitrary subset ofU with unit duration. 
There exists some constant c large enough, such that if L > c In m, then in polynomial time we can find a subset 
1Z C S and a schedule oflZ with duration at least Lj lnm, such that the remaining load Ls\R. > L/2. 

Proof. We take covers from S one by one. Let Li be the load of the remaining sensors after the i th cover has been 
taken. For the (i + l) th cover X, we take each remaining sensor into X with probability p = c In m/L,. Then we 
check if (1) A" is a valid cover, and (2) the remaining load L i+ i > Li — \ In m. For any x € U, the probability that x 
is not covered is at most (1 — p) Li < mT c , so (1) occurs with probability at least 1 — m l ~ c (probability of union of 
events). For any x G U, the probability that < Li — i lnm is at most m~( 2c_1 ) / 8c (Chernoff bound), so (2) 

occurs with probability at least 1 — TO 1 -( 2c - 1 ) 2 / 8c . Thus we can choose c large enough so that both (1) and (2) occur 
with high probability (e.g., > 1/2). We repeatedly take X until this happens and then proceed to the next cover. We 
repeat this procedure until Lj+i drops below L/2, and the lemma follows. □ 

The basic idea of our algorithms is the following. Take a partition of U with a small number of cells, and then crop 
R(-) so that each sensor fully covers a number of cells but is completely disjoint from the rest. We ensure that the load 
does not decrease by more than a constant factor and then apply Lemma 3.7. 

Theorem 3.8. If each sensor s £ 5 has unit duration, then there is a polynomial-time 0(log(n/ ' L))-approximation 
algorithm for CUBE COVER. 

Proof. We assume L > c In n for some large constant c; otherwise we just take one cover, and the theorem follows. It 
is well known that sets of rectangles in the plane admit (1/ r)-cuttings; there exists a subset TZcSofr log r rectangles 
such that in the partition A-r determined by the rectangles of 1Z, each face is intersected by the boundaries of at most 
cn/r rectangles of S [4]. We choose r = \2cn/L], so cn/r < L/2. 

Let / be a face of An, and let Sf C S denote the subset of rectangles that fully contain /. Since the load at every 
point in / is at least L and only L/2 rectangles partially cover /, we derive \S/\ > L/2. Now replace each rectangle 
R{s) by a cropped region that consists of all faces of An that s fully covers. This yields an instance of SENSOR 
COVER, with a universe of size r 2 log 2 r and load V > L/2. Applying Lemma 3.7 yields the desired result. □ 

An improved bound can be obtained when all the i?(-)' s have the same size by a more careful cropping scheme. 

Theorem 3.9. If each sensor s£5 has unit duration and each R(s) is a unit square, then there is a polynomial-time 
0(log(X max / 1 L))-approximation algorithm for CUBE COVER, where L max = max^ L(x). Note that L max < n. 
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Proof. We assume L > c In L max for some large constant c; otherwise we just take one cover, and the theorem follows. 
We draw a unit-coordinate grid T inside U. There are only 0(n/L) cells in T. For a cell 7 6 V, let £(7) C 5 denote 
the set of squares of S that intersect 7. Let n ma2: = max 7 |<S(7)|. Packing arguments imply n max < 4L max . 

Two cells in T are independent if they are at least two grid cells apart from each other in both dimensions. It is 
easy to see that we can partition T into 9 independents sets Fx, . . . , Tg, where all cells in any one set are mutually 
independent. In the following, we will show how to make fi(L/ \n(n max /L)) covers for T±, such that the remaining 
load is at least L/8. Then we repeat the process for T2, ■ ■ ■ , Tg, and ultimately we derive a schedule that covers all 
cells with duration 17 (L/ \n(n max /L)) = f2(L/ln(L max /L)). 

By the definition of independence, we can isolate the cells in Ti and only need to show that for any 7 € I\, we can 
make il(L/ \n(n max /L)) covers from 5(7) without decreasing the load of any of its neighboring 8 cells by more than 
a factor of 8. Since the load of 5(7) inside 7 is at least L, following the same approach as in the proof of Theorem 3.8, 
we can build a partition A in 7 and its neighboring cells such that each face of A is intersected by the boundaries of 
at most L/2 squares from 5(7). A has size r 2 log 2 r, where r = \2cn max / L\ . We further partition the faces of A 
that are intersected by the boundary of 7, such that each face of A is either inside 7 or outside. This increases the size 
of A by a factor at most 2. Let T be the set of faces of A that are fully covered by at least L/4 sensors from 5(7). 
T includes all faces inside 7 and some faces outside. For any face not in T, the load of S \ S(-f) must be at least 
L — L/2 — L/4 = L/4, so we can ignore it. Consider the faces in T . We crop the squares of 5(7) according to T 
in the same way as in the proof of Theorem 3.8. After cropping, by construction the load at each face of T is still at 
least L/4. Then we apply Lemma 3.7 with U = T and 5(7), which gives us f2(L/ \n{n max /L)) covers while the 
remaining sensors have load at least L/8 for any face of A. □ 

Remark. These results can be extended to any collection of shapes that admit small cuttings: disks, ellipses, etc. 



4 Sensor Cover 

Now consider the general SENSOR Cover problem, in which each R(-) is an arbitrary subset of a finite set U of size 
I U\ = 0(n). We show that a random schedule of the sensors yields an 0(log ?i)-approximation with high probability. 
This result extends that of Feige et al. [5], which deals with the unit duration case. 

Let T = cL/hxn, where c is some constant to be determined later. We show that if we choose the start time of 
each sensor randomly between and T, then we will have a valid schedule with high probability. In order to avoid 
fringe effects, we must choose positions near or T judiciously. More precisely, for a sensor s of duration d(s) < T, 
we choose its start time t(s) uniformly at random between —d(s) and T; if t(s) < 0, we reset it to 0. If d(s) > T, 
we simply set t(s) = 0. Divide T evenly into 2n time intervals [to = 0,t±], [ii,f2], ■ ■ • , [t2n~i,t2n = T], each of 
length T/2n. If d(s) > T/n, it is easy to see that for any x G R(s) and in any time interval, x is covered by s with 
probability at least (d{s) - T/2n)/{T + d(s)) > \ ■ d{s)/T. 

Consider any x € U, and let {si, . . . , s^} be the set of sensors live at x with durations at least T/n. We know that 
Yli—i d(si) > L — T/n ■ n > L/2. In any time interval [ti, ti+i], the probability that x is not covered is at most 

n(-f) s n-(-f)^(i)-(^)^-*. 

There are only 0(n 2 ) different (x, [ti, ti+i]) pairs, so the probability that some x E U is not covered at some time is 
at most 0(n 2 ) ■ n~^ = O ^n 2- ^^ . Choosing c < 1/16 yields a high probability of obtaining a valid schedule. 

The algorithm can be de-randomized using the method of conditional probability. We omit the details. 

It is not hard to see that Set Cover Packing can be reduced to Sensor Cover. Given an optimal schedule 
produced by an algorithm for SENSOR Cover, we can "snap" each starting time t(s) to the integer \t(s)] without 
introducing any gaps or decreasing the total duration. Hence, the lower bound of Feige et al. [5] applies. 

Theorem 4.1. There exists a polynomial -time 0(\ogn)-approximation algorithm for the SENSOR Cover problem. 
This bound is tight up to constant factors. 
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5 Conclusions and Open Problems 



Many questions remain open. Ideally we would like to prove stronger hardness results or find better approximation 
algorithms in order to narrow the gap between our lower and upper bounds. In fact, we have not ruled out the possibility 
of a PTAS for the Restricted Strip Cover problem, although it cannot be in terms of L. It would be interesting 
to see if other techniques for geometric optimization problems could be applied to our problem as well. 

We are also interested in understanding preemptive schedules better. For Restricted Strip Cover, a simple 
algorithm based on maximum flow yields an optimal preemptive schedule in polynomial time. In higher dimensions, 
however, it is not fully understood in which situations non-preemptive schedules are sub-optimal when compared with 
the best preemptive schedules. In general, we would like to uncover the relationship between the load of the problem 
instance, the duration of the optimal preemptive schedule, and the duration of the optimal non-preemptive schedule. 

Acknowledgement. We thank Nikhil Bansal for pointing us to the paper by Kenyon and Remila [7]. 
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A NP- Completeness of Dynamic Storage Allocation 

The following proof was given by Larry Stockmeyer to David Johnson and cited as the "private communication" 
behind the NP-completeness result in Garey and Johnson [6, Problem SR2]. To our knowledge, this proof has not 
previously appeared; we include it here essentially verbatim for the more specific results needed in Section 2.2. All 
credit goes to Larry Stockmeyer. We thank David Johnson for supplying the proof to us. 

Dynamic Storage Allocation 

Instance: Set A of items to be stored, each a E A having a size s(a), an arrival time r(a), and a departure time d(a) 
(with d(a) > r(a)), and a positive integer storage size D. 
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Question: Is there an allocation of storage for A; i.e., a function a : A — ► {1, 2, . . . , £>} such that for every a 6 A the 
allocated storage internal 1(a) = [cr(o), cr(a) + s(a) — 1] is contained in [1, D] and such that, for all a, a' € A with 
a 7^ a', if 1(a) fl J(a') is nonempty then [r(a), d(a)) fl [r(a'), d(a')) is empty? 

Theorem: Dynamic Storage Allocation is NP-complete, even when restricted to instances where s(a) e {1,2} 
for all a. 

The reduction is from the 3 -PARTITION problem. 
3 -Partition 

Instance: Set W of 3m elements, a bound B, and a positive integer size z(io) for each w G W, such that Ylwew z ( w ) — 
mB, and B/A < z(w) < B/2 for all w. 

Question: Can W be partitioned into disjoint sets W\, . . . , W m such that, for 1 < i < m, Y^weW z ( w ) = 

3-PARTITION is strongly NP-complete [6]; i.e., there is a polynomial p such that it is NP-complete when restricted 
to instances / where z(w) < p (the length of /) for all w. The condition B/4 < z(w) < B/2 is not used in the 
following reduction. 

Given an instance of 3-PARTITION as above, the corresponding instance of DYNAMIC STORAGE ALLOCATION 
has storage size D = m(B + 1) + 2. The instance is described by giving a time-ordered sequence of arrivals and 
departures of items of various sizes. It is also convenient for the description to allow an item t of size 2 to arrive several 
times, provided that it first departs before arriving again. Thus, if this item t arrives k times in the entire description, 
there are really k different items, ij for 1 < i < k, all of size 2, and no two of them exist at the same time. 

In the following, the items fi, gi, and hi all have size 1. Begin by having D items /i, . . . , /d arrive. Next fi and 
f2 depart, then t arrives and departs, then gi and gi arrive. Now do in sequence for i = 2, 3, ...,£) — 1: 

1. items gi and fi + i depart; 

2. then item t arrives and departs; 

3. then items hi and gi + i arrive. 

Finally, departs and hi arrives, and then g o departs and h o arrives. At this point, it can be seen that the order of the 
items in storage must be hi, /12, ■ ■ ■ , ho-2, hD-i, ho or hi, /12, ■ ■ ■ , hjj, ho-i or the reversal of one of those 

orderings. 

Now hi departs for every i < D — 2 that is not a multiple of B + 1. At this point, storage consists of m blocks of 
free space, each block has length B, and there are barriers (namely, hs+i, /12B+2, ■ ■ ■) between each pair of adjacent 
blocks. 

Consider first the case that sizes of items are not restricted to 1 and 2. For each w £ W, an item q w of size z(w) 
arrives. If the 3-PARTITION instance has a solution W\ , . . . , W m , then the items q w with w £ Wi can go into the first 
block of free space, the items q w with w G W2 can go into the second block of free space, and so on for W3, ■ ■ ■ , W m . 

Conversely, if all the q w fit into the m free blocks of length B, then there must be a solution to the 3 -PARTITION 
instance. 

Consider now the case that sizes of items are restricted to 1 and 2. For each w G W, the item q w of size z(w) is 
replace by z(w) items q w ,i for 1 < i < z(w), all of size 1. For each w, additional arrivals and departures are now 
added to ensure that the z(w) items q w ^ were placed in the same block of free space. The method is similar to the one 
used above with /'s, g's, and h's. First q w> i and q w ^ depart, then t arrives and departs, then q' w x and q' w 2 arrive. The 
for i = 2, 3, . . . , z(w) — 1: q' wi and q w ,i+\ depart, then t arrives and departs, then q'^ A and q' w i+1 arrive. If the q Wt i 
were not placed in the same free block, then at some point the two units of free space formed by the departure of two 
items must be separated by a barrier, so the item t of size 2 cannot be allocated. 

This is a polynomial-time transformation, because the numbers z(w) are bounded above by a polynomial in the 
length of the 3-PARTITION instance. Note further that the reduction constructs a DYNAMIC STORAGE ALLOCATION 
instance of uniform load with the ultimate determination being whether OPT = L. 
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